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We propose a relativistic gravitational theory leading to modified Newtonian dynamics, a paradigm that ex- 
plains the observed universal galactic acceleration scale and related phenomenology. We discuss phenomenolog- 
ical requirements leading to its construction and demonstrate its agreement with the observed cosmic microwave 
background and matter power spectra on linear cosmological scales. We show that its action expanded to second 
order is free of ghost instabilities and discuss its possible embedding in a more fundamental theory. 


Introduction. — Alternative theories of gravity to general 
relativity (GR) have received immense interest in the past 20 
years or so [I] [2]. The driving force behind this interest is not 
so much that gravity has not been tested in a large region of 
parameter space [3], but, more importantly, the cosmological 
systems residing in some parts of that region exhibit behavior 
from which dark matter (DM) and dark energy (DE), collec- 
tively called the dark sector, are inferred. 

While most investigations deal with DE, the hypothesis 
that the DM phenomenon is due to gravitational degrees of 
freedom (d.o.f.) has received less attention [414]. Earli- 
est evidence for the existence of DM was later sup- 
ae by observations of the motion of stars within galax- 
ies [18] [19]. Milgrom proposed [20122] that this could, 
aa result from modifying the inertia or dynamics of 
baryons or the gravitational law at accelerations smaller than 
ag ~ 1.2 x 107!°m/s?. The latter is further explored in 
where if gradients of the potential ® are smaller than ao, non- 
relativistic gravity is effectively governed by 
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Here, Gy is the Newtonian gravitational constant, and p the 
matter density. These models are referred to as modified New- 
tonian dynamics (MOND). 

Much work has gone into deducing astrophysical conse- 
quences of MOND, its consistency with data [24H46], and 
alternative DM based explanations of this law It is in- 
herently nonrelativistic and, thus, difficult to test in cosmolog- 
ical settings (but see [51]}) as systems such as the cosmic mi- 
crowave background (CMB) require a relativistic treatment. 
CMB physics involves only linearly perturbing a Friedmann- 
Lemaître-Robertson-Walker (FLRW) background, making it a 
particularly useful system, devoid of nonlinear modeling sys- 
tematics, for testing relativistic MOND (RMOND). Relativis- 
tic theories that yield MOND behavior have been proposed 
making clear predictions regarding gravitational 
lensing and cosmology. In cases where the CMB and mat- 
ter power spectra (MPS) have been computed, no theory has 
been shown to fit all of the cosmological data while preserv- 
ing MOND phenomenology in galaxies [68H76], (though see 
(77). 

We present the first RMOND theory which repro- 
duces galactic and lensing phenomenology similar to the 
Bekenstein-Sanders Tensor- Vector-Scalar (Te VeS) theory 


and, unlike TeVeS, successfully reproduces the key cos- 
mological observables: CMB and MPS. We describe its con- 
struction, discuss its cosmology and show that it is devoid of 
ghost instabilities. We discuss open questions and possibili- 
ties toward its more fundamental grounding. 

Phenomenological requirements. — RMOND theories 
have always been constructed on phenomenological grounds 
rather than based on fundamental principles. Quite likely the 
reason is that the MOND law is empirical, and even the ob- 
servation that it is scale invariant has not yet led to 
a definitive conclusion as to how this invariance could lead to 
a MOND gravitational theory. RMOND theories should obey 
the principle of general covariance and the Einstein equiva- 
lence principle. These are, however, do not provide any guid- 
ance as to how RMOND should look like. Indeed, many theo- 
ries obeying these have nothing to do with MOND, and many 
RMOND theories obeying these same principles are in con- 
flict with observations. Principle-based MOND theories in- 
clude [80H82], however, these are nonrelativistic. Still, the 
phenomenological approach, that we also follow, can provide 
valuable guidance toward a more fundamental theory. 

What are the necessary phenomenological facts that any 
successful MOND theory should lead to? It must (i) return 
to GR (hence, Newtonian gravity) when V® > ag in qua- 
sistatic situations while (ii) reproducing the MOND law (1| (ip 
when Vö < ao. It should also (iii) be in harmony with 
cosmological observations including the CMB and MPS, (iv) 
reproduce the observed gravitational lensing of isolated ob- 
jects without DM halos, and (v) propagate tensor mode grav- 
itational waves (GWs) at the speed of light. 

We consider each requirement in turn. Clearly, (i) means 
that when |V®| >> ao, the standard Poisson equation V2® = 
4rG yp holds while (ii) means that when IVd] a ao the 
MOND equation (I) holds. While in many cases [56] [60] [61] 
the transition between (i) and (ii) depends only on VO}, in 
TeVeS it is facilitated by a scalar d.o.f. y. We follow the latter 
and assume that the physics encapsulated by (i) and (ii) fits 
within the TeVeS framework. 

A template nonrelativistic action then, is 
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where ® = Ê + y is the potential that couples universally to 
matter, G is a constant and Y = |Vy|?. The field y obeys 


ior + 70] +00}, @ 


Ñ- ((dI/dy) Vy] = = 4rĜp while Ê obeys the Poisson equa- 
tion V2 = AnGp. Emergence of MOND is then ensured if 
ae wie as Vy — 0. It is in this limit that ag 
appears. 

For a point source of mass M, the MOND-to-Newton tran- 
sition occurs at ry ~ /(GnM/ao). A MOND force 
~ yGynMao/r lends its way trivially to a Newtonian force 
GnyM/r? asr < rm but in the inner Solar System this is 
not sufficient. Corrections to r~? due to y will compete with 
the post-Newtonian force ~ (Gy M)?/r?, and these are con- 
strained at Mercury’s orbit to less than ~ 1074 [83] [84]. Sup- 
pressing these may happen either through screening or track- 
ing. In the former, ọ is screened at large Voy so that ®© ~ Ê 
while in the latter y + ©/,, so that Gy = (1+1/As)Ĝ. We 
model both with A, since screening is equivalent to As — oo. 
In terms of J, tracking happens if 7 —> AsV, while screening 
occurs if J has terms Y? with p > 3/2 (this may be in conflict 
with Mercury’s orbit even as p — oo) or via higher-derivative 
terms absent from (2). 

Consider requirement (iii), that is, successful cosmology. 
In (2) we have a new d.o.f. p(T) and we expect that the same 
will appear in cosmology, albeit with a time dependence, i.e. 
o(t). Consider a flat FLRW metric so that goo = —N? and 
gij = @a?yij where N(t) is the lapse function and a(t) the 
scale factor. What should the expectation for a cosmological 
evolution of #(t) be? The MOND law for galaxies is silent re- 
garding this matter. There is, however, another empirical law 
which concerns cosmology: the existence of sizable amounts 
of energy density scaling precisely as a~°. Within the DM 
paradigm such a law is a natural consequence of particles 
obeying the collisionless Boltzmann equation. The validity 
of this law has been tested and during the time be- 
tween radiation-matter equality and recombination it is valid 
within an accuracy of ~ 1073. Do scalar field models leading 
to energy density scaling as p ~ a~? exist? 

The answer is yes: shift symmetric k essence. It has been 
shown that a scalar field with Lagrangian ~ K(4’) where 
Æ = ¢?/N?, leads to dust (i.e. p ~ a~?) plus cosmologi- 
cal constant (CC) solutions provided (4) has a minimum at 
X = Xo #0. Such a model is the low energy limit of ghost 
condensation although the latter also contains higher 
derivative terms ~ (Ld)? in its action. The FLRW action is 
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where O = ¢/N and H = a/a. Interestingly, (2) and (3) are 
shift symmetric in y and ¢ respectively. 

We propose that the MOND analog on FLRW is given by 


(3) with 
K=-2A +K2(O — Oo)? +... (4) 


where A is the CC, K and Qo parameters and (...) denote 
higher powers in this expansion. Expanding in Q — Qo rather 
than ¥ — Xo is the most general expansion leading to dust 


solutions and includes the (4) case. The CC in this model 
remains a freely specifiable parameter, just as in the A-cold 
dark matter (ACDM) model. Following [88] [89], we call this 
the (gravitational) Higgs phase. 

Requirement (iv), that is, correct gravitational lensing with- 
out DM, requires a relativistic theory. A minimal theory for 
RMOND is a scalar-tensor theory[23] with the scalar pro- 
viding for a conformal factor between two metrics. How- 
ever, since null geodesics are unaltered by conformal transfor- 
mations, such theories cannot produce enough lensing from 
baryons in the MOND regime. Sanders solved the lensing 
problem by changing the conformal into a disformal trans- 
formation using a unit-timelike vector field, incorporated 
by Bekenstein into TeVeS. The unit-timelike vector has 
component A? ~ ,/—g°9 and this ensures that the two metric 
potentials are equal (as in GR), so that solutions which mimic 
DM also produce the correct light deflection. 

Meanwhile the anisotropic scaling of the MOND law ~ 
IVg]? compared with a well-behaved cosmology implying 
terms like $2 and o, heuristically implies (gravitational) 
Lorentz violation. A good way of introducing such an ingre- 
dient is via a unit-timelike vector field A,,, much like the spirit 
of the Einstein-A&ther theory [90] 91], and TeVeS (53) (54). 

The advanced Laser Interferometer Gravitational Observa- 
tory (LIGO) and Virgo interferometers observed GWs 
from a binary neutron star merger. Combined with electro- 
magnetic observations [94], this strongly constrains the 
GW tensor mode speed to be effectively equal to that of 
light. By analyzing the tensor mode speed, TeVeS has been 
shown to be incompatible with the LIGO-Virgo ob- 
servations for any choice of parameters. The necessary d.o.f. 
@ and A, are also ingredients of TeVeS, only there, a sec- 
ond metric was introduced as a combination of gv, @ and 
A,,. In D9], ¢ and A, were combined into a timelike (but not 
unit) vector B,,, and it was shown that TeVeS may be equiva- 
lently formulated with a single metric g,,, minimally coupled 
to matter, and B,, with a noncanonical and rather complicated 
kinetic term. A general class of theories based on the pair 
{9,v, Bu} was uncovered where the tensor mode speed 
equals the speed of light in all situations, satisfying require- 
ment (v). 

The new theory. — A subset of the general class 
depends on a scalar ¢ and unit-timelike vector A“ such 


that [100] 
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where Fy, = 2V [Áv], J, = A°V.A,, and the Lagrange 
multiplier A imposes the unit-timelike constraint on A,,. In 
addition F(Y, Q) is a free function of Q = A" V,,¢ and Y = 
ge’ V ,oVid where quv = Juv + A pAr is the three-metric 
orthogonal to A”. Notice that is shift symmetric under 
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On FLRW ¢ = ¢(t) while Ao = —N and A; = 0, hence 
y = O and Q = Q. We define K(Q) = —4F(0, Q) so 
that (5) turns precisely into B}. which we have argued that it 
satisfies requirement (iii). 

In the weak-field quasistatic limit, we set goo = —1 — 24 
and gi; = (1 — 2®)yij and assume that A” aligns with the 
time direction so that A? = 1 — W and A’ = 0. The scalar is 
expanded as ¢ = 6+y with ¢ < |V¢| and ¢ may be set to its 
(late Universe) FLRW minimum Qo. Hence, Q = (1— Y) Qo. 
Then (5) leads to Y = ® which can be subbed back to get 
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where J (Y) = F(Y, Qo)/(2 — Kg). Compared with (2) a 
new term appears which looks like a “mass term” for ®, with 


H= cen Qo. The solution for ® will be as obtained from 


only for r S ro where ro ~ (ru?) ue and oscillatory 
forr = rc. We require u~! > 1Mpc so that MOND behavior 
according to may still be attained in galaxies. Thus, the 
quasistatic limit has at least three parameters: Às, ao and p. 

While matter couples only to ®, gravity comes with two 
potentials ® and p whose action is not diagonal but con- 
tains the mixing term JV ,ọ — VU - Vy. Without the 
latter, p decouples and no modification of gravity arises in 
this eee apart from u?® which is akin to ghost conden- 
sation [88] [89]. Diagonalizing by setting © = Ê + y and 
a. G = (1 — 42)G turns {6} into (2) (plus the 7? 
term). Since, UV = @, a) leads to the right lensing whenever 
the solution for ® mimics DM. This satisfies requirements (i), 
(ii) and (iv). 

Cosmological observables. — The theory just presented 
was constructed to lead to a FLRW universe resembling 
ACDM. Given a general K(Q), we define the energy den- 
sity as 8nGp = Q% — K and pressure as 8nGP = K 
so that the usual FLRW equations are satisfied. The field 
equation for ¢ may be integrated once to give 5 = iy for 


d 
initial condition Jọ. When K obeys the expansion Phen 


Q = Qo + Ip/a? + ..., so that p = o/a? +... hete 

8TGPo = Qolo. The pressure is P = wopo/a® +... where 

wo = toe is the equation of state ata = 1, that is, 
0 


w = w/a? +... so that P = wp. A time-varying w im- 
plies an adiabatic sound speed c?, = dP/dp = apse 
and if K obeys (4) then c2} = 2wo/a? +.... Clearly, w > 0 
and Cj > 0, where the zero point is reached as a —> oo. As 
the solution depends on the initial condition Jp, the density p 
is not (classically) predicted. 

For a proper cosmological matter era in the Higgs phase we 
need wy to be sufficiently small. Observations give 
w < 0.02 ata ~ 1074, hence, wọ X 2 x 10714. Mean- 
while, u~! > Mpc in order not to spoil the MOND behav- 


ior, leading to wo > ee > 1078. Unless the ef- 


fect of the jz term in (6) i is alleviated in some future theory, 
the Higgs phase cannot be extended too long in the past, and 
higher terms in must be taken into consideration. Within 
the present setup, one can arrange this with a function K(Q) 
which suppresses w and c2} during most of the cosmic evolu- 
tion. Examples are K = 2K2Z% [cosh(Z) — 1] (“Cosh func- 
tion”) and K = 2K2Zê le — 1| (“Exp function”) where 
= (Q — Qo)/Zo. 

The tight coupling of baryons to photons in the early Uni- 
verse leads to Silk damping and wipes out all small-scale 
structure in baryons, preventing the formation of galaxies in 
the late Universe. Within GR, cold DM sustains the gravita- 
tional potentials during the tight coupling period, driving the 
formation of galaxies and affecting the relative peak heights 
of the CMB as further corroborated by e.g. the Planck satel- 
lite [101]. Checking whether this theory fits the CMB and 
MPS spectra requires studying linear fluctuations on FLRW. 

We consider scalar modes in the Newtonian gauge so that 
Joo = —(1 + 2y), goi = 0 and Jij = a? (1 = 20)yi; and 
perturb the scalar as 6 = $+ y and the vector as A,, = 
{—1 — Y, Via}. The perturbed Einstein, vector and scalar 
equations, then depend on the new scalar modes y and a and 
their derivatives. The shear equation remains as in GR, as do 
the usual perturbed Boltzmann equations for baryon, photons 
and neutrinos, since they couple only to guv. 

Setting x = y + ġa, y = $ — OV, E = a+ Y and defining 
the density contrast ô and momentum divergence 0 via 


1 1 = 
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p 


0=4 
$ 
the Einstein equations take the same form as in GR, i.e. 
6G? = 8rG bar pro, and ôG’; = —8rG (Pr +Pr)V 591 
where the index J runs over all matter species including the 
new variables ô and 0. These obey standard fluid equations 
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but with nonstandard pressure contrast: 


T= 245 - —24_V? [KpE+(2-Kzs)x] (1) 


87Ga?p 
Hence, the resulting system is not equivalent to a dark fluid: 
the nonstandard pressure, thus defined, does not close under 
the fluid variables but, rather, depends on the vector field per- 
turbations a and Æ. The latter evolves with 
dk 
dQ 


+ (H + 6) Te 3caHóe] (12) 
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FIG. 1. The CMB temperature (T) C77 and E-mode polariza- 
tion CP” angular power spectra for ACDM and this theory for 
a collection of functions and parameter values. The ACDM pa- 
rameters are angular acoustic scale 1000; = 1.04171, DM den- 
sity Q.h? = 0.1202, baryon density Qh? = 0.02235, reion- 
ization optical depth 7 = 0.049, helium fraction Yye = 0.242, 
primordial scalar amplitude 10°A, = 2.078 and spectral index 
ns = 0.963, while the MOND curves deviate from these within 
~ {0.07, 0.33, 3.98, 14.29, 1.57, 0.58, 2.60} percent. MOND mod- 
els have As; = oo and their other parameters are shown in the ora i 
panel, with Qo and Zo in Mpc™t. The “Higgs-like” function param- 
eters are incompatible with a MOND limit. 


Cosmologically, the necessary additional free parameters to 
ACDM are X, (influencing the effective cosmological gravita- 
tional strength), Kg, Kə (or equivalently wo) and Qo. These 
fix u appearing in the quasistatic regime. More elaborate func- 
tions K(Q) introduce further parameters, e.g. Zo in the case 
of the “Cosh” or “Exp” functions above. Note that ag does not 
appear in the linear cosmological regime but will play a role 
once nonlinear terms from F(Y, Q) kick in. 
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where we have used the desired late Universe limit for which 
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In Figs. [1] and [2| we show the CMB and MPS in the 
case of a “Cosh”, an “Exp” and a “Higgs-like” function 
K(Q) = aay (2 - 
ing the FLRW background and linearized equations using our 
own Boltzmann code [102], which is in excellent agreement 
with other codes, see for a comparison. We have used 
adiabatic initial conditions and a standard initial power 
spectrum Py = A,k”s with amplitude A, and spectral index 
ns. The MPS has an additional bias parameter b. We used 
RECFAST version 1.5 for modeling recombination and have 
boosted k sampling, time sampling and £ sampling accuracy 


for ensuring robust results. The detailed cosmology and the 
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FIG. 2. The linear MPS P(k) for the models of Fig. |1| showing 
excellent fits to the Sloan Digital Sky Survey (SDSS) data release 7 
(DR7) luminous red galaxies (LRG) [105]. We also include a bias 
parameter b. Note that the (derived) Hubble constant for each model 
is different. 


dependence of the spectra on the parameters will be investi- 
gated elsewhere [104]. For a wide range of parameters, this 
relativistic MOND theory is consistent with the CMB mea- 
surements from Planck. This happens because c?, and w are 
small enough so that II — 0 and we get dustlike evolution as 


6 = 36 — KO and Ê = W, while the vector field decouples. 


Stability and waves. — Now, we consider stability of 
the theory on Minkowski spacetime. We expand g#” = 
nv — h”, split A, = (—1 + $h®, A;) and let Vid = 
(Qo + ġ, Vy) with h”, A; and ọ being small perturbations. 
Expanding (5) to second order gives 
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F/Q —> —2d?K /dQ? = —4K3 and OF /OQ = F = 0. 


We set OF /OY = (2 — Kpg)às as a free parameter which is 
zero in the MOND limit but nonzero in the GR limit when 
reached by tracking. Inspecting (13), the tensor mode action 
is as in GR as expected. 

For vector modes we choose the gauge h°? = 0, h°? = W* 
and h = 0 while A; = 2; and y = 0 where W* and {3° are 
transverse. Setting all modes œ e%(—“'+*-*), the dispersion 
relation for 6? is w? = k? + M? where their mass is M? = 
0 eee, hence, they are healthy if 0 < Kp < 2 and 
vs > T They decouple from T, and are not expected to be 
generated to leading order by compact objects. 

Considering scalar modes in the Newtonian gauge we set 
A = —2Ẹ, h®™ = Q and hid = —267') while A; = Via and 
find the dispersion relations w? = 0 and w? = = CK) ( + 


5K prs)k? + M?. Thus, we require that Ky > 0 in addition 
to the vector stability conditions. Only two normal modes ex- 
ist implying the presence of constraints. These are revealed 
through a Hamiltonian analysis which also shows that these 
conditions lead to a positive Hamiltonian for the 
w # 0 modes. The w = 0 case leads to a constant mode 
with zero Hamiltonian but, also, to a mode varying linearly 
with t. The Hamiltonian for the latter is positive for momenta 
larger than ~ pu and otherwise negative, also requiring that 
As > 0. Such instabilities are likely akin to Jeans-type insta- 
bilities and do not cause quantum vacuum instability at low 
momenta [108]. 

Discussion. — MOND has enjoyed success in fitting 
galactic rotation curves and reproducing 
the baryonic Tully-Fisher relation [31] [33| [45]. The radial ac- 
celeration relation (RAR) finds a comfortable interpre- 
tation within MOND. Studies of MOND with galaxy clus- 
ters report that either ag is larger in clus- 
ters and/or an additional dark component is necessary even 
when the MOND prescription is used. These studies, how- 
ever, use the classic modified-inertia MOND while the theory 
presented here has additional features warranting its separate 
testing with clusters. We note that a RAR for clusters was re- 
ported [112], similar to the galaxy one albeit with ag a factor 
of 10 higher. MOND has been tested with dwarf spheroidal 
galaxies where discrepancies for some were later dis- 
missed with improved data [28] [T13H116]. There, good agree- 
ment was reported, except for Draco and Carina where the fits 
are quite poor [26][114] [T16] [117]. It is argued that those 
two might be systems not in equilibrium. The global stabil- 
ity of M33 has been tested with positive results while 
wide-binary data do not yet yield a decisive test [119]. 

We have shown how the cosmological regime of this theory 
reproduces the CMB and MPS power spectra on linear scales 
and that MOND-like behavior emerges in the quasistatic ap- 
proximation. The latter is expected to hold for virialized ob- 
jects, however, how such objects emerge from the underly- 
ing density field, i.e. how the two regimes connect, is an 
open problem. This will happen at a scale which is expected 
to depend on ag, ys and Qo and quite likely the nonlinear 
~ V(V¢)?/ao term coming from F will play a role. It is 


reasonable to expect that on mildly nonlinear scales, the qua- 
sistatic regime is not yet reached. 

We remark that A,, also contains a pure vector mode per- 
turbation which is expected to behave similarly as in the 
Einstein-Æther theory [90] QT]. This may lead to imprints on 
the B-mode CMB polarization signal [120]. 

Pe M? 2 — 1/(87G) and canonically normalizing as 


= = SKM in (4), the FLRW action (3) becomes 
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where A? = MyV2K2 Qo. Considering the MOND 
limit in gives M2F /2 > |Vdl3 /A2 where AZ = 
12 [K2(1 + 1/As)/(2 — Kp)|*/? Mpao. This scale is indica- 
tive of the energy scale above which quantum corrections may 
be important and below which we can trust the classical the- 
ory. Since ag ~ Ho/6 then Ao = meV ~ (0. ee New- 
ton’s r7? law has been tested down to ~ 52um [121] and the 
curves in Figs|1]and[2|have Azt < 100nm. 

Absence of ghosts to oada ‘order signifies a healthy the- 
ory that could arise as a limit of a more fundamental theory. 
We do not have such a theory at present but we discuss a case 
that may bring us closer. The vector in does not seem to 
obey gauge invariance but in the quadratic action it does 
so through mixing with diffeomorphisms of h,,,,. This is not 


(14) 


an accident. Let us normalize via - — for some 


. Varying with 


A and using the constraint to foe . from the action, per- 
form a Stiickelberg transformation A, > A, + VuE/Maac 
and define the covariant derivative acting on “angular field” € 
as Da = Vué/Maac +Â. The action turns to S = SEH + 


J dta y=al — gi PH + {DMD + Mèco) } 
— WAg, 
The jee action is that of the gauged 


plus p dependent terms, where Pu = VÂ, 


M& eg., 
Teac = ~ KgM2 


ghost condensate (GGC) [122] or bumblebee field [123] 124] 
which has been proposed as a pute gauge-invariant theory 
of spontaneous Lorentz violation. The Einstein-/4ther theory, 
part of 5), is the (healthy) decoupling limit of GGC by taking 
Mecc > œ if 0 < Kg < 2 (in our notation) [122]. It is 
argued that Macc can be as high as 10!?GeV. 

Given that ¢ is shift symmetric it is natural to charge it un- 
der this symmetry similar to € letting D, = V ġ/Macc + 
Ân. Interestingly, we may identify Q— Qo > D"ED,,¢ while 
the term J’V,.¢ 4 Ê YD „ED ġ, both multiplied by appro- 
priate constants. The terms involving VY may be constructed 
using (gh” + DYED’ E/Méac) D, „oD, o. Although extend- 
ing our work as such does not explain the MOND term V3/?, 
it may provide promising directions for further improvements. 
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